Bottleneck-induced transitions in a minimal model for intracellular transport 
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We consider the influence of disorder on the non-equilibrium steady state of a minimal model 
for intracellular transport. In this model particles move unidirectionally according to the totally 
asymmetric exclusion process (TASEP) and are coupled to a bulk reservoir by Langmuir kinetics. 
Our discussion focuses on localized point defects acting as a bottleneck for the particle transport. 
Combining analytic methods and numerical simulations, we identify a rich phase behavior as a 
function of the defect strength. Our analytical approach relies on an effective mean-field theory 
obtained by splitting the lattice into two subsystems, which are effectively connected exploiting the 
local current conservation. Introducing the key concept of a carrying capacity, the maximal current 
which can flow through the bulk of the system (including the defect), we discriminate between the 
cases where the defect is irrelevant and those where it acts as a bottleneck and induces various novel 
phases (called bottleneck phases). Contrary to the simple TASEP in the presence of inhomogeneities, 
many scenarios emerge and translate into rich underlying phase-diagrams, the topological properties 
of which are discussed. 

PACS numbers: 02.50.Ey, 05.60.-k, 64.60.-i, 72.70+m 



I. INTRODUCTION 

The effects of disorder on the phase behavior has been 
investigated in a multitude of statistical mechanics mod- 
els and there is by now a good understanding of the 
ensuing equilibrium phenomena [1, 2]. In contrast, in 
non-equilibrium statistical mechanics, the effect of disor- 
der on dynamics and non-equilibrium steady state is far 
from being well understood [3] . 

One of the best studied cases is the totally asymmetric 
simple exclusion process (TASEP) [4-20]. It is defined 
as a one-dimensional lattice gas where particles are hop- 
ping stochastically in one direction subject to hard-core 
repulsion (for reviews see [21, 22]). Two types of disorder 
have been studied: hopping rates may either depend on 
the particle attempting to jump (particle-wise disorder 
[4-6] or each lattice site may be associated may be asso- 
ciated with a random quenched hopping rate (site-wise 
disorder) [7-20]. As intuitively expected, such defects 
generically induce phase separations if the particle traffic 
exceeds a certain threshold [6]. 

Our studies are motivated by an important class of 
intracellular transport processes mediated by molecular 
motors like kinesin, dyneins or myosins moving along cy- 
toskeletal filaments like microtubules or F-actin [23] . The 
dynamics of each of these molecular engines is a compli- 
cated stochastic process which we idealize by a Poisson 
process with a single rate limiting step only. There is now 
convincing evidence that transport along the cytoskeletal 
filaments is one-dimensional and binding sites are peri- 
odically spaced [24] . Since each of the binding sites can 
be occupied by at most one of the molecular motors, 
particle exclusion can play a crucial role. Though there 



is no clear evidence in vivo, at the moment, that mo- 
tor densities are large enough for particle exclusion to 
dominate the transport properties, there are in vitro in- 
vestigations underway studying the transport of kinesin 
along microtubules at high volume concentrations [25]. 
Each molecular track has a finite length, and in gen- 
eral one would like to allow for the enzyme reservoirs at 
both ends to have different densities and/or attachment 
rate at the left and right end to be different from each 
other and different from the hopping rate in the bulk. All 
this taken together defines the TASEP, introduced orig- 
inally as model for the kinetics of biopolymerisation on 
nucleic acid template [26, 27]. For it to be a proper min- 
imal model for molecular intracellular transport it also 
has to account for the fact that microtubules are embed- 
ded in the cytosol with a reservoir of motors in solution. 
This allows for motors to attach from the solution to the 
molecular track or detach from it and become part of the 
reservoir again [28] . Then one arrives at the TASEP with 
Langmuir kinetics (TASEP/LK) introduced in Ref. [29], 
which exhibits phase separation even in the absence of 
any defects. 

In this paper we would like to study the effect of an 
isolated defect on the non-equilibrium steady state of this 
minimal model for intracellular transport. We focus on a 
site-wise disorder which may be mediated by structural 
imperfections of the microtubular structure or proteins 
associated to the microtubules that change the affinity 
of the motors with the track. There has been evidence 
that these microtubule associated proteins might even be 
responsible for some diseases connected to motor proteins 
[30, 31]. 

Let us now introduce the model under investigation. 
We consider a simple exclusion process on a finite lattice 



with N sites (labeled i — 1, . . . , N), where the occupation 
number rii of each site can be either or 1. The dynamics 
of the system is described by a fully unidirectional con- 
tinuum stochastic process in which each particle jumps 
randomly with rate n to its right-neighboring vacant site. 
At the left boundary particles are introduced in the lat- 
tice with rate a, while at the right boundary they are 
extracted with rate (3. This defines what is known as 
Totally Asymmetric Simple Exclusion Process (TASEP) 
for j-j = r (we define the time scale by setting r = 1). We 
supplement this process in two ways, as shown in Fig. 1: 
(a) the system is coupled to a bulk reservoir via Langmuir 
kinetics; namely particles can attach to a site in the bulk 
with rate uja and can detach with rate ujd [32] ; (b) a spe- 
cial site k of the system represents a bottleneck and has 
a slower hopping rate r^ = q < 1 [33]. Both cases have 
been separately studied previously and show non-trivial 
features. In this work we investigate the combined effect 
of these two perturbations to the usual TASEP: our study 
concerns a detailed analysis of the non-equilibrium steady 
state properties, with emphasis on the resulting phase- 
diagrams. Two essential ingredients to carry on such an 
analysis are the local density of particles pi = (ni) and 
current profile ji = r, (n, (1 — rij+i)); where the brackets 
stand for the average over the histories. Following the 
same steps as in Refs. [29, 34] (see also [35] for a recent 
review), the stationary density is shown to obey an hier- 
archy of equations involving nearest-neighbor correlation 
functions: 

= rj_i(rij_i(l - m)) -r»(ni(l -m+i)) 

+ oo A {l-m) -u r (m), i = 2,...,N; (1) 
= a(l-ni)-(ni(l-n 2 )>; (2) 

= (n N -i(l-riN)} -/3(njv), (3) 

where, to account for the defect at site i = k, we have: 



1 if i ^ k , 

< q < 1 if i = k . 



(4) 



To be sure to capture an interesting interplay between 
boundary induced phenomena and bulk dynamics (see 
Refs. [29, 34, 36]) we shall consider a mesoscopic limit 
where local adsorption-desorption rates are rescaled such 
that the gross rates Qa,d are comparable to the injection- 
extraction rates at the boundaries: 



Hd = Nuj]j; Qa = Nuja ■ 



(5) 



Keeping fixed the gross rates, in the limit N — > oo, in- 
troduces interesting competition between boundaries and 
bulk. 

This work is organized as follows: the next section 
summarizes the main properties of the TASEP, TASEP 
coupled to the Langmuir kinetics and TASEP with a sin- 
gle inhomogeneity. In Section III, we outline the effective 
theory on which we build our analysis. Section IV is de- 
voted to the discussion of our results (phase-diagrams 
and density profiles) . Finally, we present our conclusions 
in Section V. 



t- 



UJA 
\ 



u)D 



t r 
4* 



< 1 



A 



tt 



4(- 



1 



N 



FIG. 1: Schematic representation of the TASEP with on-off 
kinetics in the presence of a bottleneck at the site i — k. The 
allowed moves are: forward jump (with rate q ^ 1 in i = k 
and r — 1 elsewhere) , entrance at the left boundary (with rate 
a) and exit at the right boundary (with rate /3), attachment 
(with rate uja), and detachment (with rate ujd) in the bulk. 



II. REVIEW OF PREVIOUS RESULTS: 

SEPARATE ROLE OF ON-OFF KINETICS AND 

DEFECT 



The effective theory presented hereafter is built on the 
properties of models akin to the one under consideration 
here, namely the TASEP and TASEP coupled to Lang- 
muir kinetics. It is thus appropriate to briefly review the 
main features of the latter and to summarize the result 
for the simple TASEP in the presence of a single inho- 
mogeneity. 



A. The TASEP 

In the absence of attachment and detachment and with 
uniform hopping rate, the model defined above reduces 
to the TASEP. Much is known about the non-equilibrium 
steady state of this paradigmatic model. Both exact 
methods [37-41] and approximated mean-field solutions 
[26, 27, 42, 43] show that as a function of the entrance 
and exit rates there are three distinct non-equilibrium 
steady states. For (3 < 1/2 and a > (3 there is a high 
density phase (HD) , where both the density and the cur- 
rent are determined by the exit rate 0. Mean-field theory 
gives a spatially constant density pi = 1 — j3 larger than 
1/2 and a constant current ji = (3(1 — (3). Thus the cur- 
rent is dominated by the low exit rate which acts as a 
bottleneck for the transport. 

In contrast, for a < 1/2 and a < (3 the low entrance 
rate is the limiting factor for the particle current which 
is now given by ji = a(l — a). Since pi — a is always 
smaller than 1/2 in this parameter range, the phase is 
also termed the low density (LD) phase. 

If both a and (3 become larger than the critical value 
1/2 the density becomes constant p* = 1/2 indepen- 
dently of the parameters at the boundaries. The cur- 
rent is limited by the particle exclusion in the bulk and 
its maximal value is j* = 1/4; therefore, this phase was 
termed maximal current (MC) phase. 



B. The TASEP with on-off kinetics 

Supplementing the TASEP, a genuine driven (non- 
equilibrium) process, with on-off (or Langmuir, equilib- 
rium) kinetics results in a system termed TASEP/LK. If 
the rates of the Langmuir kinetics are faster than, or com- 
parable to, the hopping rates, the on-off kinetics always 
dominates the driven process. To guarantee the particles 
to cover a relevant portion of the lattice before detach- 
ing, the mesoscopic limit mentioned in Eq. (5) has been 
introduced [29]. This imposes a time scale where com- 
petition between the boundary, driven processes and the 
on-off kinetics, is effective and results in rich collective 
phenomena. 

Considering Eqs. (l)-(3) with r, = 1, the mean-field 
analysis simply consists in neglecting any spatial correla- 
tions resulting in the following decoupling approximation: 

(rijnj + i) « (m) (n i+ i) = PiPi+i ; (6) 

ji ~ K)(l - (n i+ i)) = Pi(l -Pi+i)- (7) 

Taking the continuum limit with the new spatial variable 
< x = i/N < 1 and using the mesoscopic limit of 
Eq. (5) one obtains: 



(2p - l)d xP - Q DP + Q A {1 - p) = . 



(8) 



For the sake of simplicity, let us consider the case where 
the two rates are the same (£Id = ^A = fi). In this case, 
Eq. (8) reads: 



(2p-l)[d xP -n}=Q. 



(9) 



This equation obviously has two solutions: a linear one 
with slope fl and a constant one coinciding with the crit- 
ical density of the TASEP: 



p(x) = fix + C and p(x) = p* 



(10) 



where C is a constant to be determined by the boundary 
conditions. The linear density profile results in a space- 
dependent current, which in mean-field reads 



j[p{x)} = p(x) [1 - p(x)} 



(11) 



Equation (9) has to be supplemented by the boundary 
conditions 



p(0) = a; p(l) = l-(3. 



(12) 



When the solution of Eq. (10) cannot be matched contin- 
uously with the left and right boundaries (12), the den- 
sity profile displays a localized discontinuity (or shock) 
in the bulk. This translates into the emergence of mixed 
phases. The latter are discussed in Refs. [29, 34] and are 
summarized in the TASEP/LK phase-diagram of Fig. 2a. 
Let us illustrate this concept by considering the transi- 
tion from the LD to the LD-HD phase. In the former, 
the density is determined by the left boundary and reads 



Lowering the exit rate (3 (with fixed a), there is a site 
where both the currents imposed by the left and the right 
boundaries, meet at x w (see Fig. 2c). This is a mixed LD 
and HD phase. As shown in Fig. 2b, the density has a 
sharp phase boundary (shock) at x w : 



p{x) 




Qx + a 

1 - /? + fl(x - 1) 



for < x < x v 
for x w < x < 1 



(14) 



Another feature of the TASEP/LK is the particle-hole 
symmetry, i.e. the properties of the system are invariant 
under the exchanges a <-> j3, x <-> 1 — x, and p = 1 — p. 
This translate in Fig. 2 which is symmetric w. r. t. the 
linea = /3 [29]. 

The case Qa ^ &d follows along the same lines but 
is mathematically more involved; for its treatment we 
refer to Ref. [34] . In the following we will consider Ha = 
Qd = &> and present the case Ha ^ &d only at the end 
of Section IV. 
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FIG. 2: (a) Phase diagram of TASEP/LK for K = 1. One rec- 
ognizes seven phases: in addition to the TASEP LD, HD and 
MC phases; there are four more coexistence phases, namely 
the LD/HD, LD/MC, MC/HD and LD/MC/HD phases. The 
shaded region highlights the LD/HD coexistence where a lo- 
calized domain wall appears, (b) Typical density profile in 
the LD/HD phase and (c) the corresponding current profile. 
At the matching point x v between the left (j a ) and right (jp) 
currents a domain wall develops and connects the left (p a ) 
and right (pg) density profile. 



p(x) = p a (x) — fix + a . 



(13) 



The TASEP with a single inhomogeneity: a 
brief review 



Before presenting the details of our theoretical treat- 
ment of the TASEP/LK system in the presence of a bot- 
tleneck, and to gain some intuitive understanding of the 
underlying physics, it is convenient to outline the prop- 
erties of the simple TASEP perturbed by a localized in- 
homogeneity [7, 10]. Here we consider the TASEP with 
a defect at site k where the hopping rate is q < 1. 



Consider the LD phase, where the system is 'diluted' 
and the particles are well separated. In such a "low traf- 
fic" situation, one does not expect any macroscopic ef- 
fects arising from the presence of the bottleneck. One 
rather expects a local peak in the density profile (see 
Fig. 3a): pi^k — ct and pk = a + h (LD phase), where 
h is the height of the local jump imposed by the defect. 
Since an exclusion process without coupling to a bulk 
reservoir has a spatially constant current ji = a(l — a), 
one finds h = h(a, q) = a(l — q)/q- Thus the height of 
the density peak increases with a and q^ 1 . This can 
certainly not happen without bound. In fact, h can- 
not exceed ft max = 1 — 2a. This may be seen as fol- 
lows. If ft, > /i max , there would be a site i\ such that 
p h < a + /i max = 1 - a and p il+1 > a + h max . This 
case has to be discarded as it contradicts current conser- 
vation: ji 1 = pi t (1 — pi x ) < ji = a(l — a). Therefore, for 
given a and (3, there is a critical value of g* below which 
the local peak is no longer a possible solution. The crit- 
ical g* follows from the requirement that the height of 
the peak cannot exceed ft max , i.e. h(a,q* a ) — h max . One 
thus finds 



9 Q = 



1 



(15) 



Equivalently, if one keeps q fixed, a local peak is obtained 
for 



a < 



1 + q 



(16) 



It directly follows from the underlying particle-hole 
symmetry that the same reasoning holds in the HD phase 
(a > (3 and (3 < 1/2). There the system is so 'packed' 
that the bottleneck is only responsible for a local dip in 
the density profile: pi^k = 1 — (3 and pk = 1 — (3— h{(3, q), 
with h(/3, q%) < h' max = 1 - 2/3. Thus, in the HD phase, 
a local dip is only possible for 



9/3 



or, keeping q fixed, for 



3 



1-/3 



<q 



1 + 9 



= 0c 



(17) 



(18) 



We conclude that, for a fixed value of q, there are criti- 
cal values for the entrance and exit rates above which the 
defect leads to a jump in the density profile. In contrast 
to the local density perturbation ("spike") for small a 
and /3, this is a macroscopic effect and does not vanish in 
the thermodynamics limit (while the spike width scales 
as 1/JV). 

Denoting respectively by p L and p R the densities on the 
left and right of the inhomogeneity, one expects a sharp 
step profile displaying a jump at the defect position: 



S and p ' 



(19) 
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FIG. 3: Transition between a spike (a) and a step (b) in 
the density profile for the TASEP with defect. Stochastic 
simulations are performed on systems of size N = 128 with 
parameters q — 0.1, Xd = 0.5, (a) (a,/3) = (0.05,0.8), (b) 
(a, /3) = (0.8,0.8). 



Thus, the current through the defect is 



jk = qpk(i -p/b+i) = q 



6 



1 



, (20) 



while the current flowing through the bulk of the track 

reads 



Jijtk = 



1 



1 



-6 



(21) 



As the current is conserved, these expressions should co- 
incide, which gives the jump in the density profile: 



25 



1-g 
1+9 



(22) 



This is an important quantity which measures the 
strength of the defect: S — when q = 1 and 5=1/2 for 
3 = 0. 

Thus, for the density profile, one finds 

P L = -^— and P R = ^t, (23) 

1 + q q+ 1 

while the q— dependent current flowing through the en- 
tire system, playing the role of an effective (bottleneck 
induced) MC, reads 



3d = 



(1 + q) 2 4 



(24) 



Since there is a jump and the density profile is flat on 
both sides of the defect, one may effectively split the sys- 
tems into two parts connected by current conservation 
[10]; (see Fig. 5). The effective exit rates of the left sub- 
system and the entrance of the right one are therefore: 



C^cff — /3eff — 



1 



(25) 



It follows from this discussion that for the TASEP only 
the MC phase is affected by the presence of a bottle- 
neck, while the LD and HD phases remain unaltered. 
The resulting phase-diagram therefore still displays the 
same topological features (with three phases: LD, HD 
and MC) as in the homogeneous case, with the lines 
a = (3 = 1/2, delimiting the MC phase, lowered to by 
a = a c and /3 = j3 c . 



III. EFFECTIVE MEAN-FIELD THEORY 

In this section, we describe an effective mean-field 
theory for the non-equilibrium steady state of the 
TASEP/LK in the presence of a bottleneck. 
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FIG. 4: Transition between spike (a) and step (b) for 
the TASEP/LK with defect. Stochastic simulations are per- 
formed on systems of size N — 1 28 with parameters q = 0.1, 
x d = 0.5, n = 0.1, (a) («,/3) = (0.05,0.8), (b) (a,fJ) = 
(0.8,0.8). 

As for the TASEP perturbed by the presence of a lo- 
calized defect, there are regions of the parameter space 
where the inhomogeneity affects the system only locally. 
In other regions, the defect has macroscopic effects and 
is said to be relevant (see Fig. 3). In this case, the 
properties of the system are studied by splitting the 
lattice at the defect site k into two subsystems (again 
termed L and R) and performing a continuum limit (see 
Fig. 5). In such a limit, the position of the defect becomes 
x d = limjv— >oo k/N, and the density can be written as 



_ j p L {x) ,0<x<x d 
P[X) \p*{x) ,x d <x<l 



(26) 
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FIG. 5: Schematic representation of the division into two sub- 
systems allowing to apply a mean-field theory. The last step 
(bottom of the figure), illustrates the continuum limit that 
we are considering (see the main text). 

In the presence of a bottleneck, the current is locally 
conserved for the same reasons as in the TASEP/LK sys- 
tem [34] [essentially because the attachment/detachment 
rates scale as in Eq. (5)]: 



3i+\ - 3i 



A^ 



(Qa + Qd)p 



i+lj 



JV- 



->0. (27) 



This allows to couple both subsystems R and L along 
the same lines as for the simple TASEP perturbed by 
an isolated inhomogeneity, with effective right and left 
boundaries given by Eq. (25). 

Suppose now that we are in the parameter range where 
the defect dominates the current and density profile. 
Then we may argue as in the previous section and ob- 
tain for the current at the defect site : 



3k = Jfc+i = 



(1 + 9) 2 
and the corresponding densities 



1 



Pk 



i + q 



Pk+i = 



l 



(28) 



(29) 



Again this implies a defect induced jump in the density 
of magnitude 



28 



1 



Pk - Pk+i 



l + « 



(30) 



A key difference to our previous discussion of a pure 
TASEP in Sec. II is that these properties apply only lo- 
cally in the vicinity of the defect. In the continuum limit, 
the density profile is linear [see Eq. (10)] and has to match 
the boundary conditions given by Eqs. (29). Hence the 
density profile p d {x) imposed by the defect reads 



Pd{x) 



'si( x -x d ) + ^- (L) 

n(x- Xd ) + T f- (R) 



(31) 



Since the mean-field current-density relationship is given 
by Eq. (11), this immediately implies for the current 



jd = Q\x - x d \ 



<l 



l + « 



1 



1 



— filiEd — X 



(32) 



In stark contrast to the simple TASEP, here the current 
is a space-dependent quantity. As illustrated in Fig. 6, 
the signature of a defect is the depletion on a macroscopic 
scale of the current profile (Fig. 6b). Correspondingly, 
the density displays a 'zigzag' shape with a jump at x d 
and a linear profile in its vicinity as expressed in Eq. (31) 
(see Fig. 6a). 

At x — x d ± £, the current imposed by the defect 
reaches the maximal value j* — 1/4. Thus, the deple- 
tion or screening length £ induced by the bottleneck is 
the solution of pd(x d ± £) = p* = 1/2 and reads 



*-i 



(33) 



The defect is thus screened best for a strong coupling to 
the reservoir (fi 3> 8) and, of course, if the inhomogeneity 
is weak, i.e. <5 — ► 0. 

Since the screening length increases with the strength 
of the defect it can even become larger than the length of 
the sub-systems. This happens when £ is larger than the 
distance to either one of the boundaries. According to the 
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FIG. 6: Sketch of the signature of the strong defect (zigzag 
shape) in the density (a) and current (b) profile. 



conditions Xd — £ < and Xd + £ > 1, this happens when 
5 > max((Si,^2) (and fixed 0) or, equivalently, when 17 < 
min(£li, f^) (with fixed q), where 



(5i = flxd 
_S_ 

Xd 



fii 



and 
and 



S 2 = 0(1 - Xd ) (34) 

0, ' 



1 - Xd 



(35) 



If <5 < miri(<$i, #2) the screening length is shorter than 
both lengths of the subsystems. When the location of 
the defect is not centered, two additional cases arise: for 
<5i < S < &2 (#2 < 6 < 6i) the screening length stays 
within the subsystem R (L), while it is larger than the 
size of the sub-lattice L (R). 

In the TASEP and in the TASEP/LK the maximal 
current that can flow through the system is a constant 
j* = 1/4. For the TASEP with a defect that value 
is lowered to j* d = q/(l + q) 2 [sec Eq. (24)]. As the 
current is space dependent and locally conserved in the 
TASEP/LK with an inhomogeneity, the maximal flow of 
particles through the bulk varies spatially. This suggests 
to term such a quantity the carrying capacity C(x) of the 
system. 

The more drastic effect of the bottleneck appears when 
its strength satisfies S > max(<5i, 62)- In this case, defect 
screens the entire system and, as shown in Fig. 7a, the 
carrying capacity reads 



C(x) =Ci(x) = jd{x). 



(36) 



In intermediate cases, when the screening length cov- 
ers part of the system, one has three possible scenarios 
illustrated in Fig. 7b- 7d, namely 



C(x) = C 2 (x) 



C(x) =C 3 (x) 




, \x-Xd\ <£ 
, else 

< x < Xd + £ 
x d + £ < x < 1 



(37) 



(38) 



C(x) = Ci(x) = 



(39) 



j* , < x < x d - £ 

[jd(x) ,x d -£ < x < 1 

For each possible carrying capacity of Fig. 7, one ob- 
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FIG. 7: The four typical carrying capacity profiles C(x) 
displayed by the system (parameters Q = {0.3,1.5,0.5,0.5}, 
Xd = {0.5,0.5,0.7,0.3}, q = 0.3). Depending on the defect 
strength q, on the position Xd of the defect and on the rate Q, 
the defect imposed current jd combines in four different ways 
with the maximal current j*. Each of these profiles induce 
topologically distinct phase diagram (see the text). 

tains different phase diagrams. All these scenarios are 
discussed in the next section. 

To determine the global current and density profiles 
from the carrying capacity, one compares the latter with 
the currents imposed by the open boundaries: 

j a = (a + ilx) {1 - a - ilx) , (40) 

Jp = \p + Sl(l-x)][l-0-il(l-x)], (41) 

corresponding to the left and right density profiles 



p a (x) = fix + a, 

PP {x) = n(x-i) + (1-/3). 



(42) 
(43) 



As C(x) acts as an effective maximal current in the bulk, 
j a (x) and jp{x) cannot exceed its value. Actually, as 
for the TASEP (where C — 1/4), the entrance and exit 
boundaries only matter on the macroscopic part of the 
system where they impose currents smaller than C(x). 
Thus, by matching j a (x) and jp(x) with the carrying ca- 
pacity, one determines boundaries separating boundary- 
induced phases and mixed ones. 



IV. RESULTS: PHASE-DIAGRAM OF THE 

TASEP/LK IN THE PRESENCE OF A 

BOTTLENECK 

In this section we discuss the four possible scenarios 
arising for each carrying capacity presented above. We 



explicitly construct the density profiles and the phase dia- 
grams when the screening length is larger than the size of 
the two subsystems (carrying capacity C\{x)) and when 
£ is shorter than the size of the two subsystems (carry- 
ing capacity C2(x)). In the last subsection we extend the 
results obtained so far for the case of equal attachment 
and detachment rates (Qa = ^d) to the more general 
situation £Ia ^ 91d- 



A. Large screening length: the case C\{x) 

Let us first consider the case where the carrying capac- 
ity is entirely determined by the defect, i.e. C\ (x) = jd(x). 
In this situation, sketched in Fig. 7a, the bottleneck is 
strong enough and always imposes a current jd{%) < j* ■ 
Then we may distinguish between three cases depending 
on the magnitude of the current imposed by the left and 
right boundary (see Fig. 8). 
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FIG. 8: Construction of the current profile. Top: For 
given carrying capacity C\{x) (solid line) and fixed exit rate 
j3 — 0.2, various scenarios arise when the entrance rate 
a — {0.05, 0.2, 0.6} is varied. Bottom: Three different pro- 
files are shown in the small graphs (solid line: global current 
profile, dashed line: defect and boundary currents) emerging 
from three different left boundary conditions. Parameters are 
x d = 1/2, q = 0.3, fi = 0.3, (3 = 0.2, a = {0.05, 0.2, 0.6}. 

There are two extreme cases. The entrance and exit 
boundary currents j a and jp exceed the carrying capac- 
ity C\{x), the current through the system settled at the 
value jd(x) imposed by the defect (this is situation for the 
left subsystem presented as case 3 in Fig. 8). Thus, the 
density exhibits the piecewise profile given by Eq. (31). 
As the current is independent of the left and right bound- 
aries, this phase is termed pure bottleneck phase (BP). 



In contrast, for low entrance and exit rates one recovers 
the TASEP/LK density profile perturbed by a local spike 
or a dip (see previous section). This is presented as case 
1 in Fig. 8. 

In intermediate regions of the parameter space, a sit- 
uation like the one presented as case 2 in Fig. 8, ap- 
pears. Similarly to what happens for the homogeneous 
TASEP/LK, when the densities p a , pp and pd cannot be 
matched continuously, shocks form in the density profile 
and then we have coexistence of several phases. This 
can happen either on the left or right subsystem. The 
positions of the shocks follow from the local conserva- 
tion of the current, i.e. j a (x^) — jd(%w) ( on L) and 
i/s(a£) = idOO (° n R )- With Eq 8 - (32)-(41), one finds 



Xd 



1 
9. 



2\ 1 + Xd ■ n 



(3 + 6 



(44) 
(45) 



The conditions for having domain walls within the two 
subsystems are < x\, < x d and x d < x^ < 1. This 
translates into the following conditions for the entrance 
and exit rates: 



1 



5 — Qx c ] < a < 

2 2 

--6-n(l-x d )< (3 < - 



5 + Qx d , (46) 

5 + n(l-x d ). (47) 



These conditions allow to identify two values of the 
boundary rates a~ and /3~ for which a domain wall en- 
ters respectively the left or right subsystem and two val- 
ues a£ and /3+ for which the domain wall enters the left 
or right subsystem: 



1 



a c = - - 5 ± flx d , 
Pt = ±-6±fl(l-x d ). 



(48) 
(49) 



Within the range defined by these critical rates, a shock 
is localized in each subsystem. At values above u~ and 
j3~ the defect becomes relevant and at least part of the 
density and current profile is dictated by the defect. This 
suggests to term the corresponding regions of the param- 
eters space bottleneck phases. 

Let us now construct the density profiles in the var- 
ious bottleneck phases focusing first on the left subsys- 
tem. Depending on the entrance rate, one distinguishes 
two cases: (i) The density profile corresponding to case 
2 in Fig. 8, arising when a ~ < a < a+, is presented in 
Fig. 9a (numerical simulations are discussed later). Here, 
j a intersects the left branch of C\(x) at x]^, where a do- 
main wall forms. Then the density (of the left subsystem) 
reads: 



L , . J p a (x) ,0<x<x. 



p d (x) , x^ < x < x d 



(50) 



This coexistence phase is called LD-BP since it is char- 
acterized by the coexistence of a low density and a bot- 
tleneck phases, (ii) The density profile corresponding 
to case 3 in Fig. 8, arising for a > a+, is shown in 
Fig. 9b. Here, j a is always above the left branch of 
C\{x) and the defect imposes an effective high density 
phase (called simply BP) corresponding to the maximal 
current jd on subsystem L. The corresponding density 
profile reads p L (x) — pd(x). 

One proceeds in a similar way for the right subsys- 
tem (R). When /3~ </?</?+, there is coexistence be- 
tween high-density (boundary induced) and an effective 
low-density (defect induced) phases, called BP-HD (see 
Fig. 9d). This results in a domain wall at x|J and in the 
density 



Rf n _ I Pd{x) , x d < x < x\ 

\Pp{x) . ^w < X < 1 



(51) 



S > max{<5i, 62}- We have checked the predictions of 
our MF theory against stochastic numerical simulations 
(following the Bortz-Kalos-Lebowitz scheme for kinetic 
Monte Carlo [44]) and, as shown in Fig. 9, have found 
good agreement (both qualitative and quantitative) with 
the predictions of the MF theory. Of course, due to finite- 
size effects, boundary layers form in the vicinity of the 
shocks and 'soften' the transition in the density profile. 
When the system-size is increased, the accuracy of the 
MF theory improves and the density profile displays a 
sharp jump. In addition to boundary layers, as already 
noted in Ref. [7] for the simple TASEP perturbed by one 
defect, in the bottleneck phases, one remarks a slight 
but systematic deviation from the MF predictions (see 
e.g. Figs. 9b and 9d). This effect is due to correlations 
not taken into account by the mean-field approximation 
and is long-range, scaling as the inverse of the distance 
to the defect [7] . 



When j3 > /3"J~, the defect dominates and imposes an 
effective low density phase, called BP where the density 
is p R (x) = pd(x). 

Considering the whole system, by combining the above 
mixed phases (on R and L) we obtain the following four 
bottleneck phases: 
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FIG. 10: Phase diagram for Q = 0.3, q = 0.4 and x d = 1/2, 
i.e. for the carrying capacity C\(x) (large screening length). 
Continuous lines are the phase boundaries introduced by the 
defect; dashed lines are the phase boundaries already present 
in the model without bottleneck. The shadowed region in- 
dicates the bottleneck phases where the defect is relevant, 
the darkest one highlights the pure bottleneck phase (for the 
meaning of the different phases see text). 



FIG. 9: Examples of density profiles in the bottleneck 
phases for a carrying capacity C(x) = Ci(x). Stochastic sim- 
ulations (continuous line) are compared to analytical mean- 
field predictions (dashed line). The system size is N = 4096 
and the parameters are q = 0.3, x d = 0.5, fio = 0.3, LD- 
BP: (a, /3) = (0.3,0.6), BP: {a, 13) = (0.6,0.6), LD-BP-HD: 
[a, (3) = (0.3,0.3) and BP-HD: (a,/3) = (0.6,0.3). 

The emergence of these four defect dominated mixed 
phases is the most dramatic effect of the bottleneck when 



The various phases are summarized in the a//3-phasc- 
diagram of Fig. 10. In this figure, because Xd = 1/2 and 
Qa = &d = Q, the particle-hole symmetry (which still 
holds in the presence of the bottleneck) results in the 
invariance of the phase-diagram w. r. t. the line a — (3. 

As can be seen from the phase-diagram of Fig. 10, 
the bottleneck-induced mixed phases (LD-BP-HD, LD- 
BP, BP-HD and BP) occupy the upper right part of 
the diagram (shadowed region in Fig. 10). Only at the 



borders of the phase-diagram, corresponding to particu- 
larly low/high entrance/exit rates, one recovers the same 
phases as in the defect-free TASEP/LK model (there the 
defect is irrelevant). For a > a+ and (3 > /3+, i.e. in the 
right top corner of the phase-diagram (darkest shadowed 
region in Fig. 10), the entire system is in a pure bottleneck 
phase. By tuning the strength of the defect <5 and the on- 
off parameter $7, the phase boundaries can be shifted to 
recover the short screening length case (discussed in the 
following) and eventually the usual TASEP/LK behavior 
(when the defect is irrelevant). On the other hand, by 
increasing the strength of the defect (or by reducing the 
on-off rate) the phase boundaries move toward the axes 
of the phase-diagram and can even merge with them. 
Actually, this occurs for a defect strength 6. Above this 
threshold, the whole a//3-phase-diagram is characterized 
by bottleneck phases (the homogeneous LD, HD and LD- 
HD phases are squeezed out of the diagram) . The critical 
strength S is therefore the maximum between the value 
determined by the conditions a~ = and (3~ — 0, which 
leads to 



Eq. (53), corresponding to a density profile 



.' 1 x 1 

d = max < Oi , — 

2 '2 



(52) 



Note that for a too strong defect, 5 > S, the control of 
the system from the boundaries is lost. 



B. Short screening length 

We now consider the case where the screening length 
is short so that the carrying capacity reaches the value 
j* = 1/4 in the bulk of the system. 



1. Symmetric screening 

We first consider the case illustrated in Fig. 7b, the 
carrying capacity loses the trace of the defect in the sys- 
tem, at a distance larger than the screening length: 



C 2 (x) 



jd{x) 



,\X- 

, else 



x d\ < S, 



(53) 



This situation arises when 5 < min{<5i, ^2} (with fixed f2) 
or, equivalcntly, when Q > max{£7i,£72} (with fixed S). 
This carrying capacity corresponds to the richest case in 
terms of new bottleneck phases, due to the profile of C2 
characterized by four distinct regions. 

As in the previous situation, the bottleneck is rele- 
vant and induces new mixed phases when a > a~ and 
[3 > (3~ , where these critical values are again given by 
Eqs. (48). Elsewhere in the parameter space, the homo- 
geneous TASEP/LK profiles locally perturbed by a spike 
(or a dip) are recovered. 

When the system is driven above its carrying capac- 
ity (a,/3 > 1/2), it exhibits the current profile given in 



f(x) = 



Pd{x) 



,Xd 
,Xd 



x > £ 
x < £ 



(54) 



As the current profile of Eq. (53) reaches the MC value, 
contrary to the case of a large screening length, this is 
no longer a pure bottleneck phase but corresponds to a 
MC-BP-MC phase (see dark shadowed region in Fig. 12). 

Between these two extremal situations, more intricate 
mixed phases appear. As in the previous case, phases 
characterized by shocks appear when the boundary cur- 
rents match jd- 

Contrary to the previous situation, where jd covered 
the whole subsystem, now the defect imposed current 
extends up to a distance equal to the screening length 
from the bottleneck. Hence, shocks can emerge only in 
a macroscopic region, Xd — £ < x}^ < Xd and Xd < x R < 
Xd + £, in the vicinity of the defect. These conditions 
translate into new critical values of the entrance and exit 
rates, namely 



1 



1 



— flXA 



Q 



P'c = 



1 



1 



Q 



- n{i - x d ) (55) 



As a new scenario, here the boundary currents can 
reach the upper-bound j* = 1/4 on the subsystems L 
and R when j a {x L ) = jfj(x R ) = 1/4. With Eqs. (40) and 
(41), one finds that this occurs at 



l-2a 



1 



1-2/3 

2n 



(56) 



We note that these quantities are independent of the 
properties of the defect, which is thus screened. Let us 
consider a density profile which exemplifies the richness 
of this case. When a' c < a < a* = 1/2 the resulting 
current j a {x) saturates at the value 1/4, while the den- 
sity reads p h (x) = 1/2 for x L < x < Xd — £. Similarly, 
on the subsystem R, when f3' c < (3 < (3* = 1/2 and the 
current jp(x) saturates at the value 1/4, the density is 
p h (x) — 1/2 for Xd + £ < x < x R . Within the range of 
the screening length, Xd — £ < x < Xd + £, the carrying 
capacity and the current flowing through the system co- 
incide with jd(x) and the density is given by Eq. (31). 
Summarizing, in this case the density profile is piecewise 
and one distinguishes five regions: 



p(x) = 



' p a (x) 


, < x < x L 


p* 


, x L < x < Xd — £, 


Pdix) 


, x d - £ < x < x d + S, 


P* 


, Xd + £, < x < x R 


.Pp{x) 


, x R < x < 1 



(57) 



This case, denoted LD-MC-BP-MC-HD, is illustrated in 
Fig. 11 and corresponds to the coexistence of a low- 
density and high-density (on L and R, respectively), 
two maximal current (one in both subsystems) and a 
bottleneck-induced mixed phases. 
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FIG. 11: Examples of density profile (top) and current (bot- 
tom) for the system in the LD-MC-BP-MC-HD phase: the 
defect is relevant and the carrying capacity is C2(x). stochas- 
tic simulations (continuous line) are compared to analyti- 
cal mean- field predictions (dashed line). The system size is 
N = 4096, K = 1, ft = ft D = 1.5, q = 0.3, (a, 13) = (0.2.0.2) 
and Xd = 1/2: one can clearly distinguish the various phases 
and note a discontinuity in the density profile in the proximity 
of the defect, while the current profile is continuous. 



The three possible scenarios for the bottleneck phases 
on the two sub-lattices results in nine mixed phases on 
the whole system, as summarized in Tab. I. 

We have also checked our MF predictions against 
stochastic numerical simulations, as illustrated in Fig. 11. 
Again, we have found qualitative and quantitative agree- 
ment. The small deviations from the MF theory observed 
in Fig. 11 can be explained along the above discussion on 
the role of the correlations and finite-size effects. 

The a//3-phase-diagram corresponding to a system 
with a carrying capacity C2(x) is shown in Fig. 12 and 
characterized by the transitions lines corresponding to 
the critical values a~,a! c ,a* and (3~ , (3' c , /?* . The nine 
bottleneck-induced phases appear in the core of the phase 
diagram (shadowed region in Fig. 12). Contrary to the 
phase-diagram obtained with C = Ci, instead of a pure 
bottleneck phase, here one finds a MC-BP-MC coexis- 
tence phase, which is independent of the entrance and 
exit, boundaries (right top corner of Fig. 12). This 
phase, obtained when a > a* and (3 > (3* , is the 'bot- 
tleneck analogous' of the MC phase in the homogeneous 
TASEP/LK system. 
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FIG. 12: Phase-diagram for ft = 0.3, x d = 1/2 and q = 0.8, 
i.e. for the carrying capacity €2(3;) (short screening length). 
Continuous lines are the phase boundaries introduced by the 
defect; dashed lines are the phase boundaries already present 
in the model without bottleneck. The shadowed region indi- 
cates the bottleneck phases where the defect is relevant, the 
darkest one highlights the bottleneck phase independent of 
the entrance (a) and exit (/?) rates (see the text and Tab. I). 



2. Asymmetric (partial) screening 

To conclude the study of the case SIa = &d, let us 
briefly consider the scenarios where the system displays 
a carrying capacity of type C3 (a;) or C4 (x) . This is possi- 
ble when the defect is not at the center of the system. In 
fact, the carrying capacity of the system is C^{x) when 
5\ < S < 82 and is C±{x) when 82 < 5 < 5\. As one 
can infer from the profile of C(x) — C^{x) [see Fig. c], 
on the subsystem L (R) this case is identical to that dis- 
cussed for C(x) — C\(x) (C(x) — C2(xj). Therefore, it 
directly follows from the above tables that there are six 
new phases for the system displaying a carrying capacity 
Cs(x). The latter are summarized in the following table: 
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Similarly, six bottleneck phases are also obtained for 
a carrying capacity C{x) = C&{x), as it is clear from the 
table below: 
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Left|/Right-> 


BP-HD 


BP-MC-HD 


BP-MC 


LD-BP 


LD-BP-HD' 


LD-BP-MC-HD 4 


LD-BP-MC 1 


LD-MC-BP 


LD-MC-BP-HD S 


LD-MC-BP-MC-HD 5 
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MC-BP-HD a 


MC-BP-MC-HD 


MC-BP-MC d 



TABLE I: Bottleneck phases for the case C(x) = £2(2;) 
part) of the phase-diagram represented in Fig. 12. 



For each phase, the label (1-9) refers to a given region of the (shadowed 



The corresponding phase-diagrams directly follow from 
those discussed above for the carrying capacities Ci(x) 
and C2(x). 



C. Topological features 

The above discussion has shown that the sole presence 
of a localized bottleneck is responsible for the emergence 
of new (sub-)phases and drastic topological changes of 
the a//3-phase-diagram. Here, we aim to discuss further 
the important structural changes induced by the bottle- 
neck in the density and current profiles of the TASEP/LK 
by considering the a/<5-phase-diagram (see Fig. 13). 

We have already seen that the most appealing prop- 
erties of the phase-diagrams of the model under consid- 
eration are the new bottleneck (sub)phases. As sum- 
marized in Figs. 10 and 12 a large portion of the a/ (3- 
phase-diagrams is dominated by the defect properties. 
The phase boundaries separating the usual phases from 
the bottleneck (i.e. a = a~ and j3 — j3~) are straight 
lines since they depend only on the detachment rate and 
on the strength and position of the defect, but not on the 
entrance exit rates. 

In the simple TASEP, the transitions between the usual 
phases and the bottleneck phases can be considered as 
discontinuous, since the influence of the defect changes 
abruptly from a local peak to a global step (see Fig. 3). 
In the case of TASEP/LK with bottleneck, a shock en- 
ters continuously the system and the transition can be 
consider indeed continuous, as it was the in the simple 
TASEP/LK (dashed line in Fig. 13). 

The phase boundaries a+ (/3+), for S > 5i (S > q^), 
and a* = 1/2 (/?* = 1/2), for S < Si (6 < 5 2 ), separate 
the point where the density profile is dominated by the 
defect from the one where the boundary currents still 
play a role. Considering the position of the domain walls 
x^ and x^ as the order parameters, the transitions are 
all continuous. A similar transition is also identified by 
the boundary a' c (/3' c ) since the matching point (no longer 
a shock) x L (x R ) enters the system continuously. 

The physical modifications in the system following the 
transition between a carrying capacity C\{x) and C^ix) 
can be observed if we plot the cut of the phase-diagram 
in S and a (choosing j3 > 1/2), as done in Fig. 13. 
In this graph, the line a~ identifies the transition be- 
tween relevant an irrelevant defect, while cu+, a' c and 
a* = 1/2 the transitions between bottleneck phases. The 
line a = 1/2 — rip was present in the pure TASEP/LK. 
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FIG. 13: Cut of the phase-diagram in <5 and a (parameters: 
Q. = 0.4, Xd — 1/2, j3 > 1/2). Solid lines identify usual phase 
transitions, while the dashed line the peculiar transition be- 
tween relevant and irrelevant defect. The multiple coexistence 
points V and Q are shown in the graph (see main text). 



The line 5 = 8\ identifies the continuous transition be- 
tween the two different carrying capacities, C\{x) and 
C2{x). There arc two relevant points in the phase di- 
agram where different phases coexist, indicated with V 
(five phases) and Q (four phases) in Fig. 13. 



D. Results and phase-diagrams when Qa 7^ f2o 

In this Subsection we shortly discuss the case where 
the attachment and detachment rates differ, Qa ^ &d- 
This is mathematically more tedious, as for the homoge- 
neous TASEP/LK model, at MF level, one needs to solve 
Eqs. (8,12) which solution implies multivalued (Lam- 
bert) functions with two real branches [34]. Introducing 
the binding constant K = Qa/Qd 7^ 1, the Langmuir 
isotherm pi reads pi = K/{K + 1). Here, having split 
the problem in two subsystems one has to consider the 
equation 



(2p 



L,R 



l)d xP ^\x) = (K + 1)Q D [ PI - p L ' H (x)] , (58) 
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again supplemented with the (subsystems) boundary 
conditions Eqs. (12) and (25) 



p L (0) = a; p L (x d ) = 



1 



p R M = 



9+1 



1 + q 
p R (l) = l-/3- 



(59) 
(60) 



We refer the readers to Ref. [34] for a detailed mathemat- 
ical treatment of this kind of equations, and report our 
results for the phase-diagram and density profiles of the 
TASEP/LK model in the presence of a bottleneck when 
fin y^ Ha (i-c K ^ 1). As for the homogeneous model, 
one can take advantage of the underlying particle-hole 
symmetry to restrict the discussion to the case K > 1 
[34] . Except for the mathematical treatment of the MF 
bulk equation, we follow the same lines as in the case 
K = 1. Again, one has to distinguish the case where the 
carrying capacity coincides with the current imposed by 
the defect (C(x) = jd(x)) from the situation where C(x) 
reaches the maximal current value j* K . While for K — 1 
the maximal current available in the bulk is j* = 1/4, 
here j* K {x) < j* is a non-constant space-dependent quan- 
tity. When C(x) — jd(x), except some topological asym- 
metries, one essentially recovers the same phase-diagram 
as in the K = 1 situation when the carrying capacity is 
of C\(x) type: as illustrated in Fig. 14, the new bottle- 
neck (sub-)phases are BP, MC-BP, LD-P and LD-BP-HD 
(see Fig. 15 (left) for an example of current and density 
profiles). Again, we notice that the core of the phase- 
diagram is entirely determined by the defect (shadowed 
region in Fig. 14), which is also responsible for a pure BP 
for sufficiently high rates a and (3. 

On the other hand, when K ^ 1 it follows from the 
analytical solution of the above-mentioned MF equation 
that there is only one screening length on the subsystem 
R. In fact, it turns out that the left branch of jd{x) can 
never reach the maximal current j* K {x) on the subsys- 
tem L(so, for if/ 1, there is no finite screening length 
on the subsystem L): the screening length may only be 
shorter than the length of the subsystem R, with the 
point Xd + £ < 1. It follows that carrying capacities of 
types Cz(x) and C±(x) are topologically prohibited (even 
in the asymmetric case x& ^ 1/2) when K ^ 1. 

In contrast to the case K — 1, as in general the cur- 
rent jp(x) can not reach the value j* K {x), the carrying 
capacity of type Cs(x) gives rise to four bottleneck phases 
(instead of six as for the case K = 1), namely the LD- 
BP-HD, LD-BP-MC, BP-HD and BP-MC. An example 
of BP-MC current and density profiles is shown in Fig. 15 
(right). 



V. CONCLUSION 

This work has been devoted to the study of the influ- 
ence of a bottleneck (point-wise disorder) on the station- 
ary properties of a biologically inspired stochastic trans- 
port model obtained by coupling two paradigmatic equi- 
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FIG. 14: Phase diagram for K = 2, Q. D = 0.1, x d = 1/2 and 
q — 0.3, i.e. for the Ci(a;)-like carrying capacity. Continuous 
lines are the phase boundaries introduced by the defect (BP); 
dashed lines are the phase boundaries already present in the 
model without bottleneck. The shadowed region indicates the 
bottleneck phases where the defect is relevant, the darkest one 
highlights the pure bottleneck phase. The numbers stem for 
the different phases: LD-BP (1), BP (2), LD-BP-HD (3) and 
BP-HD (4). For the description of the phases see the text and 
Figs. 15. 



librium and non-equilibrium processes. Namely, we have 
considered the competition between the totally asymmet- 
ric exclusion process (TASEP) and Langmuir kinetics 
(LK) in the presence of open boundaries and a bottle- 
neck, which locally slows down any incoming particles. 
The current and density profiles in the non-equilibrium 
steady state have been investigated analytically via an 
effective mean-field theory built on splitting the lattice 
into two subsystems. Our analytical results were checked 
against numerical (Monte-Carlo) simulations. 

As a consequence of the competition between the 
TASEP and LK dynamics, the effects of a single bot- 
tleneck in the TASEP/LK model are much more dra- 
matic than in the simple TASEP [10], or the so-called 
^-TASEP (TASEP for extended objects) [14], where a 
localized defect was shown to merely shift some tran- 
sitions line in the phase-diagram, but do not affect its 
topology. Here, new and mixed phases induced by the 
bottleneck have been obtained. 

As a key concept of our analysis, we have introduced 
the carrying capacity, which is defined as the maximal 
current that can flow through the bulk of the system. 
In contrast to the simple TASEP the spatial dependence 
of the current, caused by the Langmuir kinetics, makes 
the carrying capacity non-trivial: The defect depletes the 
current profile within a distance that we called screen- 
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FIG. 15: Two examples of density profile and current in the 
LD-BP-HD for the Ci(:r)-like carrying density capacity (left) 
and BP-MC phase for the Cz{x) (right). Stochastic simula- 
tions (continuous line) are compared to analytical mean-field 
predictions (dashed line). The system size is N — 4096 and 
the parameters are q — 0.2, x d — 1/2, K — 2, flo — 0.3, 
(a, /3) = (0.1.0.2) (left); q = 0.3, x d = 1/2, K = 2, Q D = 0.3, 
(a,/3) = (0.8,0.8) (right). 



ports the maximal current, independently of the bound- 
aries. Between these two extreme situations, we have 
found several coexistence phases, where the density pro- 
file exhibits stable shocks and kinks. Indeed, above some 
specific parameter values the phase-diagram is charac- 
terized by bottleneck phases. Depending on the screening 
length imposed by the defect, which can cover the entire 
system or part of it, different phases-diagrams arise. The 
latter are characterized by four, six or nine bottleneck 
phases, which have been quantitatively studied within 
our mean-field theory. 

The analysis carried out in this work can be straight- 
forwardly extended to many variants of the TASEP/LK 
model. As an example, let us mention the case of a lat- 
tice gas where dimers (modelling the usual two heads of 
molecular motors) would move as bound entities accord- 
ing to the ^-TASEP (with 1 = 2) [45] and could expe- 
rience an Langmuir-like on-off kinetics. Recently, such 
a system has been studied (without disorder) within an 
appropriate mean-field-likc scheme [36]. The point- wise 
version of this system could be investigated along the 
same lines described in the present work. 

In addition to direct extensions, we think that the 
method outlined in this paper could pave the way to 
study the TASEP/LK models in more 'realistic' and bio- 
physically relevant situations, as in the presence of clus- 
ters of competing defects or quenched site- wise random- 
ness. 



ing length. This quantity increases with the strength of 
the defect and decreases with the attachment-detachment 
rates. The competition between the current imposed at 
the boundaries and the one limited by the defect deter- 
mines the density profiles and the ensuing phase diagram. 
When the boundary currents dominant, the phase be- 
havior of the defect-free system is recovered. Also, above 
some critical entrance and exit rates, the system trans- 
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